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A detailed analysis to find a missing i/'CdS) is made by utilizing the recent precise measurements of the cross 
section for the process e*e^ —» il/(2S)n*Tt^ by Belle. Assuming three resonances y(4360), K(4660), and i^(4S) 
to fit the data, we obtain the resonance parameters for as m = 4243 MeV and T = 16 ± 31 MeV showing 
a narrow state as predicted before. A combined fit to the data —> i//(2S)7t*n^, h^n^n^, and^^ow is also 
performed to obtain the similar resonance parameters of i/'(45). The upper limit of the branching ratio is fitted 
to be S(if/(4S) —> if{2S)n'it^) < 3 x 10“^, which can be understood by hadronic loop contributions within a 
reasonable range of parameters. In addition, the ratios of the branching ratios of the >f{4S) dipion transition 
to that of i/'(45) —> are fitted, which can be further measured by BESIII and the forthcoming Bellell to 

confirm the existence of iJ/{4S). 

PACS numbers: 14.40.pq, 13.66.Bc 


I. INTRODUCTION 

In 2014, the Belle Collaboration released new experimen¬ 
tal data of the cross section of e^e~ fi{2S)n^n^ [1], where 
two charmonium-like states 7(4360) and 7(4660) were con¬ 
firmed again. Besides 7(4360) and 7(4660), the experimental 
data of the cross section for 4i{2S)n^7t^ [1] show the 

existence of an extra narrow structure near 4.2 GeV, which 
inspires our interest in exploring its origin. 

Since this peak is just located at the invariant mass of 
fi{2S)n'^7r~ from the annihilation [2], we can conclude 
that it must have the j’^^ = 1 quantum number. This fact 
leads us to an important issue to be checked, whether this 
structure around 4.2 GeV in the invariant mass of ^(2S)n^7t^ 
has a relation with the proposal of a missing xjjiAS) in our two 
recent papers [3, 4]. 

In Ref. [3], we once predicted a missing higher charmo- 
nium i/'(45) with the mass 4263 MeV when applying the sim¬ 
ilarity between the 7/i/r and T families, which is supported by 
former theoretical calculation of the charmonium spectrum by 
considering the screened potential [5-7]. Further study of the 
open-charm decay channels of i/^(45 ) indicates that \j/{4S ) has 
very narrow decay width [3], which is the reason why iA(45) is 
still missing in the present experimental data. BESIII reported 
the measurement of the cross section for —> jt^n^hc at 

^fs - 3.90 ~ 4.42 GeV [8]. Yuan performed a fit of the avail¬ 
able experimental data of —> n^n^hc from 3.90 to 4.42 

GeV, and found a narrow structure around 4.2 GeV, where the 
mass and width are reported to be M = 4216 + 7 MeV and 
r = 39 +17 MeV or M = 4230+10 MeV and F = 12+36 MeV 
[9], depending on the different assumptions of the line shape 
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trend above 4.42 GeV. The authors of Ref. [3] attributed this 
narrow structure existing in jt^n~hc to the predicted 

missing i/t(45). 

Later, the present authors noticed new results of e^e~ 
XcQ^^ given by BESIII [10], and indicated that this observed 
resonance structure with mass M - 4230 + 8 MeV and width 
F = 38+12 MeV in Xc 0 <^ [10] can be explained as 

the predicted \p{4S) [4] since the calculated branching ratio 
of ^{4S) —» uXcQ can overlap with the experimental data in 
a reasonable parameter range. In addition, the upper limit of 
a branching ratio of xjjiAS) —» ^Iso predicted to be 

1.9 X 10“^, which can be further tested by BESIII, Belle and 
the forthcoming Bellell'. 

Besides e^e~ jt^n^hc [8] and [10], 

we have been trying to search for any signal of the pre¬ 
dicted i/'(45) in other hidden-charm decay channels from the 
e^e- annihilation. The new experimental data of 
il/{2S)jt^ [2] provides us good opportunity. In this work, we 
first fit the observed cross section for e^e~ —> >f/{2S )7r^7r“ by 
including three resonance contributions, i.e., those of the pre¬ 
dicted tA^(45), 7(4360) and 7(4660). In this way, we extract 
the resonance parameters of these three resonances. Com¬ 
paring the obtained resonance parameters of if{4S) with the 
theoretical results, we can further test the \p{4S) assignment 
to this structure around 4.2 GeV in the invariant mass spec¬ 
trum of i//{2S)n'^7T^. After fitting the experimental data of 
—> il/{2S)Tt^n^ [2], we next obtain information of the 
branching ratio of f/{4S) —» \l/{2S)n^Tt'~. In this work, we 
also calculate the branching ratios of ijj{4S) ^{2S)n^7t^ 
and make a comparison with the extracted data, which can 
provide an extra test of the proposal in the present work. The 
detailed calculation will be given in the next section. 


' In Ref. [11], the branching ratio limit B{iJ{4S) —» T]Jlf) < 1.3% was 
extracted by analyzing the experimental data, which is consistent with our 
theoretical prediction in Ref. [4]. 
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Since enough experimental data have been accumulated for 
the hidden-charm decay channels from the annihilation, 
i.e., e^e~ [8], e^e" XcO(^ [10]^ ^nd 

4i{2S)n^n^ [2], we are able to perform a combined analysis 
of all these data assuming the existence of i/'(45) around 4.2 
GeV. This is also the main issue in this work. 

This paper is organized as follows. In the following sec¬ 
tion, a detailed analysis of fitting the —> i/'(25 )7r^7r“ data 

is presented on the premise that the predicted i/'(45), T(4360) 
and 7(4660) are present. A theoretical estimate of the branch¬ 
ing ratio for i/'(45) —» ^{2S)n^n^ is presented in Sec. III. In 
Sec. IV, a combined fit to the hidden charm decay channels is 
performed and Sec. V is devoted to a short summary. 


II. FIT OF e*e- if/(2S DATA 


When looking at the cross section for ^(2S)n^n^ 

[1], one can find a number of events near 4.2 GeV other than 
the structures of 7(4360) and 7(4660). Here, we perform a fit 
to the cross section for —> il/{2S)7r^7T^ with a coherent 

sum of three resonances, which are 7(4360), 7(4660) and a 
resonance around 4.2 GeV named 7(4230) in this work. The 
total cross section can be depicted by 


cr(m) = 



PS2->3M 

PS2-^3('«/) 


( 1 ) 


where is the phase angle between different resonances with 
(po - 0, and PS 2^3 indicates the phase space of 2 —> 3 body 
process. The indices / = 0,1,2 are assigned to the resonances 
7(4230), 7(4360) and 7(4660), respectively. The concrete 
form of the Breit-Wigner function of a resonance with mass 
iriR and width T^ is 


Ji 27 Tr;^-m ^ /)r« 

BW(m) = ^ -. (2) 

+ iniKi r 


TABLE I: The parameters determined by fitting the experimental 
data of e*e^ —» if/(2S)n*n^. The symbol / in the branching ratio 
indicates the final state, ili{2S )n* tt^ . The masses and total decay 
widths are in units of MeV, while the product of branching ratios 
and the dileption decay width S(R —> /) is in units of eV. 



Sol. I Sol. II Sol. Ill Sol. IV 

^^(4360) 

4356 ± 8 

ry(4360) 

65± 10 

rn:360)®(y(4360) ^ /) 

5.8 ± 0.8 7.5 ± 1.5 6.6 ± 2.5 8.6 ± 1.2 

'”y(4660) 

4656 ± 24 

ry(4660) 

73 ±29 

ryW66o,®(y(4660) ^ /) 

1.5 ± 0.5 1.6 ± 0.6 5.1 ± 0.7 5.2 ± 0.7 

</>i 

0.5 ± 0.5 0.8 ± 0.5 4.7 ± 0.2 5.0 ± 0.2 

^^(4230) 

4243 ± 7 

ry(4230) 

16±31 

P,;^230)S(7(4230) ^ /) 

0.5 ±0.5 1.4 ±0.5 0.5 ±0.4 1.4 ±0.5 

02 

5.5 ±1.0 1.7 ±0.5 5.4 ±1.0 1.6 ±0.5 

A^ndf 

45.3/64 



(GeV) 


yi (GeV) 


In the Breit-Wigner function, resonance parameters and a 
product of a branching ratio T^ !B{R ^ \lj{2S)n^n^) are 
treated as free parameters. 

Since we include three resonances to fit the cross section 
with 11 free parameters, there are four independent solutions 
with equally good quality. The obtained parameters are listed 
in Table I. With these parameters, ;t'^/ndf is estimated to be 
45.3/64. The corresponding total cross sections and reso¬ 
nance contributions are presented in Fig. 1. The masses 
and widths for 7(4360) and 7(4660) obtained from our fit 
are ?riy(^ 4360 ) — (4356 ± 8) MeV, 77 ^ 4350 ) — (65 ±10) MeV, 
nry( 4660 ) — (4656 + 24) MeV and ry( 466 o) ~ (73 + 29) MeV, 
respectively. The fitted resonance parameters of 7(4360) and 
7(4660) are consistent with those found by PDG [12]. As for 
the resonance 7(4230), the mass and width are fitted to be, 

»Jy( 4230 ) = 4243 ± 7 MeV, 

ry( 4230 ) = 16 ± 31 MeV. 


FIG. 1: (color online). Our fits of different solutions to the cross 
section data for 7(4230) —> if/(2S)n*n^. The red solid, cyan dashed 
curves and blue dots with error bars are the total results, resonance 
contributions and experimental data, respectively. 


Here, we notice that the resonance parameters of the new 
resonance 7(4230) are consistent with those observed in the 
cross sections for n^n^hc and coXcO [9, 10]. 

In addition, in our previous work [3], we have predicted a 
narrow charmonium state xpi^S) close to 4.2 GeV. The res¬ 
onance parameters of the structure in the cross sections for 
7T^7r^ip{2S) by the Belle Collaboration, 
n^n-hc {m = 4230 ± 10 MeV, T = 12 ± 36 MeV or 
m = 4216 ± 7 MeV, T = 39 ± 17 MeV) and ^ ojXcO 
(m = 4230±8 MeV, T = 38± 12 MeV) by the BESIII Collab¬ 
oration are in line with our expectation of the missing i//(4-S). 
From these fact, we suspect that this new resonance 7(4230) 
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FIG. 2: (color online). A comparison of the fits to the cross sections 
for e*e^ —> with different schemes. 


around 4.2 GeV in the cross sections of e^e~ Ji^n^ip{2S ) is 
a good candidate for our predicted i/'(45). In the next section, 
we will continue to test this point. 


It should be noticed that in the vicinity of 4.2 GeV, there 
are three resonances, i/'(4160), F(4260), and the introduced 
y(4230) in the present work. To test the the significance 
of these three resonances in the cross sections for —> 

7r^7r"i/f(25), we fit the cross sections in four different schemes; 
(a) y(4360) + F(4660), (b) 7(4260) + 7(4360) + 7(4660), (c) 
(A(4160) + 7(4360) + 7(4660), and (d) 7(4230) + 7(4360) + 
7(4660). The resonance parameters of the ^/r(4160), 7(4260), 
7(4360) and 7(4660) are constrained by the PDG values. A 
comparison of the fitting curves in different schemes is pre¬ 
sented in Fig. 2. The dominant difference of these schemes 
comes from the description of the data around 4.2 GeV, where 
the fitting curves are presented in the subfigures of Fig. 2. 
The for the different schemes are 51.2, 48.8, 49.2, and 
45.3, respectively, which indicates that the experimental data 
in the vicinity of 4.2 GeV of the cross sections for —> 

favor a narrow structure 7(4230). 


Besides the resonance parameters of 7(4230), we also ob¬ 
tain the product of the branching ratio to e^e^ifr(2S) and the 
partial width (abbreviated as FS), which can provide us 
with further more information about this state and is presented 
in Table I. Since 7(4230) is a narrow state and far away from 
7(4660), interference between 7(4230) and 7(4660) is very 
weak and ignorable, while interference between 7(4230) and 
7(4360) is significant. FS can be divided into two groups ac¬ 
cording to interference between 7(4230) and 7(4360). For 
Sols. I and III, this interference is constructive. The center 
value of FS in these two solutions is 0.5 eV, which is compa¬ 
rable to zero when including the errors. For Sols. II and IV, 
this interference is destructive, in which case FS is (1.4 + 0.5) 
eV. From our fit, we obtain the upper limit of FS to be 1.9 eV. 



FIG. 3: (color online). Typical meson loop contributions to if/(4S) —> 
if/(2S )n^n^, where the dipion comes from cr meson. 


III. MESON LOOP CONTRIBUTIONS TO 

i/r(4S) ^ >lj{2S)n*n- 

From our fit to the experimental data of e^e~ —> 
^{2S)ji^n^, we find that the product of the branching ratio to 
n^n^4t{2S) and the partial width, F^4230)®(^(4230) —> 
il/{2S)ji^n^) (FS), is dependent on interference between 
7(4230) and 7(4360). If interference is constructive, FS is 
comparable to zero, whereas if interference is destructive, FS 
is 1.4 + 0.5 eV, which indicates that the upper limit of FS is 
1.9 eV. As we have discussed in the previous section, the reso¬ 
nance parameters of 7(4230) are consistent with the expecta¬ 
tions of the missing i/'(45). To further check the possibility of 
7(4230) as the (&(45), we need to understand the fitted prod¬ 
uct of the branching ratio to n^n^ip{2S ) and the e^e~ partial 
width, which is less than 1.9 eV. Because the measurement of 
the i/'(45) dilepton decay width is not available at the moment, 
we adopt the theoretical estimate of which is 0.63 keV 
[6] or 0.66 keV [7] in the quark potential model. 

In Ref. [6], authors interpreted the observed i/r(4160) as a 
^/'(45) state and the ^^^(4415) as a state to solve the con¬ 

tradiction observed in dilepton decay widths between experi¬ 
mental and theoretical data. However, with this assignment, 
the ifr(2D) state would not have any experimental counterpart. 
We notice that the theoretical mass of il/i4-S) was estimated 
to be 4274 MeV [6], which is about 120 MeV higher than 
the mass of ^!f(4160) and is comparable with the 7(4230) dis¬ 
cussed in this work. Thus, the t/r(45) predicted by the screened 
potential model in Ref. [6] is more likely to be 7(4230) exist¬ 
ing in —> n^7T^il/{2S) rather than i/r(4160). 

With the theoretical value of we can roughly es¬ 

timate the upper limit of the branching ratio for (//{A-S) —> 
ij/{2S)n^7i^ to be 3.0 x 10“^ with the assumption of 7(4230) 
as the missing \j/{AS ). The interpretation on this branching ra¬ 
tio in the (A(45) framework can further test the nature of this 
state. 

Here, one should specify that the mass of the missing (A(45) 
is above the threshold of a pair of charmed mesons, which 
means the xjjiAS) dominantly decays into a charmed meson 
pair [3]. A charmed meson pair can transit into a charmo- 
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nium and a light meson in the hnal state, this mechanism, 
named the meson loop mechanism, plays a cracial role in 
understanding the hidden charm/bottom decay behaviors of 
a higher heavy quarkonium [4, 13-17]. The process i/'(45) 
decaying into il/{2S)n^n^ occurs via a charmed meson loop 
as shown in Fig. 3. Here, we adopt the effective Lagrangian 
approach to estimate meson loop contributions to the decay 
i/^(45) —» 1/^(25)71^71“. The effective interaction between char- 
monium and a charmed meson pair can be constructed in the 
heavy quark limit [18-20]. The specihc effective Lagrangians 
for are 

dp -Ddp D*J) 

-Dl (4) 

where Tf*'* - The effective coupling 

constants between tl/(2S) and the charmed meson pair can 
be related to the decay constant of i/'(25) by g^( 2 s)DD - 
g^{2S)D'D'mD>hno - g^(2s)D'Dmf ylmolmD^ - m^{2S)lU(2S)- 
The i/f(25) decay constant can be evaluated from the lep- 

tonic decay width of ip{2S). With F^^s) = 2.36 keV, we have 
fijjiis) - 297 MeV. As for the interaction between and 

charmed meson pairs, we adopt the same formula as shown 
above for il/{2S), but the related coupling constants are evalu¬ 
ated by the corresponding decay widths. Since the experimen¬ 
tal measurements for the open charm decay of iJ/{4-S ) are not 
available yet, here we adopt the theoretical predictions based 
on the quark pair creation model [3]. As shown in Ref. [3], 
the open charm decay width depends on the parameter R in¬ 
troduced in the spatial wave function of il/{4-S). Thus, the cou¬ 
pling constants estimated from the partial decay width will 
depend on the parameter R. The values of the coupling con¬ 
stants related to i/'(45) have been given in our previous work 

[4]. 

The effective Lagrangians related to the cr meson are [21- 
24], 




£ 


(jnn 


-ga-OO^dl^cr + ga-0'D-£>*£>*^ CT, 


(5) 


where g^-oD = go-D-D' = mo^gn! V6 with g„ = 3.73 [25, 26]. 
The coupling of crnn is evaluated by the decay width of the 
cr meson, where, Fj°‘ ^ 1'^- In the present 

work, we take nia- - 526 MeV and Fo- = 302 MeV [27]. 

In the triangle diagram, the exchanged charmed mesons are 
off-shell. To describe the off shell effect and the structure of 
the exchange mesons and to regularize the divergence, we in¬ 
troduce a form factor in the amplitudes. In the present work, 
we adopt a monopole form form factor as [4, 13, 14], 


T{q\K^) 


— A^ 
^2 - A2 ’ 


( 6 ) 


where q and are the momentum and the mass of the 
exchanged meson, respectively. The parameter A can be 


reparametrized as A = niE + a^J^QCD with Aqco = 220 MeV. 
The dimensionless parameter aA is of order 1 and dependent 
on the specihc process [28]. 

With the above preparations, we parametrize the amplitude 
of ^/r(45)(po) ^ ip(2S)(px)n^{p2)7T^{p2) in the form. 


^ ~ ^( 4 s£^{ 2 S)(fs Sfiv + foPlfiPOv) 

^_ Sa-JTTT _ 

iP 2 + p-iY -ml. + 


(7) 


where fs and fo are evaluated from the loop integral of the 
triangle diagrams listed in Fig. 3. For the broad resonance cr, 
we introduce a momentum-dependent decay width as [27] 


— Ter 


TWq- 

Tyijiji 


\p{m„„)\ 

\pima.)\ 


( 8 ) 


where \p{m„„)\ - is the pion momentum in the 

mother particle rest frame. With the above formula, we get 
the differential partial decay width as 

where the overline above |A4|^ indicates the average over the 
spin of initial i/7(45). We can estimate the partial decay width 
by integrating over m^(^ 2 S)K and m^„. 

The meson loop contributions to the branching ratio for 
(//{A-S) —> i//(2S)7r^7T^ depend on both parameters R and a a, 
which are introduced in the wave function of i//{4-S ) and the 
form factors in the amplitudes. The R and a a dependence 
of the branching ratio for i/'(45) ^ ip{2S)n^n^ is presented 
in Fig. 4. Some contour lines with several typical values of 
the branching ratio are also shown. We hnd that the branch¬ 
ing ratio S(i/f(45) ^ il/{2S resulting from meson loop 

contributions overlaps with the upper limit, 3.0 x 10“2, ob¬ 
tained by htting the cross section for ^ i/f(25)7r^7r“. 
This fact further supports the assignment of T(4230) existing 
in the ^ i/f(25 )7r^7r“ process as the missing (A(45). 


IV. COMBINED FIT TO e+e“ ^ iA(25)7r+7r“, /2,7r+7r“, XmOJ 

As we have indicated in Section II, the resonance param¬ 
eters of the new structure, T(4230), obtained by htting the 
experimental data of ^ ip{2S are consistent with 
those obtained by htting the cross sections of —> hen^n^ 

and e^e~ XcO^i which indicates that the structures in these 
three processes may come from the same source. To further 
test this conjecture, we shall perform a combined ht to the ex¬ 
perimental data of the cross sections for —> i/f(25)7r^7r“, 

—> hcTT^TT^ and ^ Xc 0 <^ simultaneously. 

We adopt the same formula as Eq. (1) to describe the cross 
section for —> Xc 0 <^ by replacing the 2^3 phase 

space with a 2 ^ 2 phase space. As for ^ 
we notice that the experimental data for the cross section of 
—» hcTT^n^ above 4.4 GeV are not available at present 
[9]. Here, we adopt two schemes to describe the cross section 
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R (GeV)"' 


FIG. 4: (color online). The R and ttA dependence of the branching 
ratio for tf/(4S) —> if/(2S)7T*n^. 


for —> hcn^n^ depending on the tendency above 4.4 GeV 

[9]. In the hrst scheme, the cross section for —> hcU^n^ 

goes down above 4.4 GeV, where two Breit-Wigner functions 
are adopted to depict the cross section as Eq. (1). In the sec¬ 
ond scheme, the cross section goes up above 4.4 GeV, where 
the cross section is described by the coherent sum of a Breit- 
Wigner function and the 2^3 phase space in the form. 


cr„ 


■^hcTi'^n ~ 


c VPS 2 ^ 3 (»t) + e''^BW(m). 


PS2^3(»t) 


I PS2_3(»*«) 


( 10 ) 


A. Scheme I 



yS (GcV) ^ (GeV) 


FIG. 6: (color online). The different solutions of the resonance con¬ 
tributions and our fitting results for the cross section for —> 

in Scheme I. The cyan dashed and red solid curves are the 
resonance contributions and the fitting results, respectively. 



FIG. 7: (color online). The different solutions of the resonance con¬ 
tributions and our fitting results for the cross section for 
(solid curve) in Scheme I. The dashed curve is the phase space of 
Xcqoj- 



•501 . II . 

4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.44.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 


^ (GeV) x/? (GeV) 

FIG. 5: (color online). The different solutions of the resonance con¬ 
tributions and our fitting results for the cross section for —> 

in Scheme I. The cyan dashed and red solid curves are 
the resonance contributions and the fitting results, respectively. 


In this scheme, we simulate the cross section for e^e~ 

and Acow by three res¬ 

onances, two resonances and one resonance, respectively. A 
resonance near 4.2 GeV, F(4230), is included in all three pro¬ 
cesses. The parameters determined by fitting the experimental 
data are listed in Table II. With these parameters , the;\f^/ndf 
is estimated to be 52.2/81. The resonance parameters of the 
7(4230) included in all three processes are fitted to be 


»Jy(4230) =4234 + 5 MeV, 
rr( 4230 ) = 29 + 14 MeV. 

The resonance parameters of 7(4230) obtained by a combined 
fit are consistent with those determined by fitting the cross 
sections for —> hcn^n^ and e^e~ [9. 10] sepa¬ 

rately. In addition, the narrow width obtained is also consis¬ 
tent with our expectation of i/'(45) in Ref. [3]. 

Besides 7(4230), two other resonances 7(4360) and 
7(4660) are also involved in the cross section for 
il/{2S)ji*n~. The resonance parameters of these two 
charmonium-like states are determined to be niYi^^ao) - 4359+ 
7 MeV, 7^(4360) = 64 + 11 MeV, my( 4660 ) = 4666 + 
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TABLE II: The parameters determined by fitting the experimental data of e*e^ —> il/(2S )n*n^, h^n^n^, simultaneously, where the 

experimental data of are depicted by two Breit-Wigner structures. The masses and the total decay widths are in units of MeV, 

while the product of the branching ratios are in units of eV. 


final State 

l//(2S )7T*K 

h^n^n 

XcOOJ 

Sol. A 

Sol. B Sol. C Sol. D 

Sol. 1 Sol. 2 


™7(4230) 

4234 ± 5 



77 ( 4230 ) 

29+14 



7^230)^^45)^/) 1.3 ±0.5 

0.3 ±0.2 1.3 ±0.5 0.3 ±0.3 

0.2 ±0.1 7.1 ±2.9 

2.2 ±0.6 

™7(4300) 


4294 ± 11 


77 ( 4300 ) 


201 ± 55 


r7(4MO,®(5'(4300) ^ /) ••• 


14.7 ±2.0 23.9 ±2.4 


Pi 


5.7 ±0.8 3.7 ±0.1 


^ 7 ( 4360 ) 

4359 ± 7 



77(4360) 

64+11 



^7^4360)^(^(4360) ^ /) 7.4 ± 1.4 5.5 ± 1.9 8.9 ± 1.0 6.6 ± 1.0 



/2 4.2 ± 0.4 

1.5 ±0.9 4.4 ±0.4 1.7 ±0.6 



^ 7 ( 4660 ) 

4666 ± 28 



77 ( 4660 ) 

90 ±20 



rn4M)®(3'(4660) ^ /) 1.9 ± 0.8 

1.8 ±0.7 6.0 ±3.2 5.8 ±2.3 



03 5.2 ± 0.7 

2.2 ± 1.0 3.1 ±0.5 0.1 ±2.1 



T^ndf 

52.2/81 




TABLE III: Our estimate of the ratios of the branching ratios of the y(4230) dipion transitions to the one of T(4230) —» XcO^- 



fntp(2S)K'^n 



Sol. A Sol. B Sol. C Sol. D 

Sol. 1 S 0 I .2 

Scheme I 

0.56 ±0.25 0.13 ±0.12 0.59 ± 0.26 0.14 ±0.13 

0.07 ±0.05 3.21 ±0.97 

Scheme II 

0.06 ± 0.09 0.26 ± 0.30 0.27 ± 0.20 0.05 ± 0.08 

0.11 ±0.04 


28 MeV, r 7 ( 4660 ) = 90 + 20 MeV, respectively, which are 
consistent with the corresponding PDG average values [12]. 
When simulating the cross section for \p{2S)n^n^ 

by three Breit-Wigner functions, we get four different so¬ 
lutions with equally good fit quality, whose parameters for 
different solutions are listed in Table II. Among these four 
solutions, the masses and widths of y(4230), 7(4360) and 
7(4660) are the same, but the product x S«^i/,( 2 S);r+;r- 
with R - (7(4230), 7(4360), 7(4660)) and the phase angles 
(pi are different. For Sol. A and Sol. C, interference between 
7(4230) and 7(4360) is destructive and FS = 1.3 + 0.5 for 
7(4230), while for Sol. B and Sol. D, interference is con¬ 
structive. In Fig. 5, the resonance contributions from 7(4230), 
7(4360) and 7(4660) are presented for different solutions. For 
Sol. A and Sol. C, the contribution from 7(4230) is signifi¬ 
cant, while for Sol. B and Sol. D, it is relatively small. 

As for the cross section for hcTT^n^, it is simu¬ 

lated by two Breit-Wigner functions for 7(4230) and 7(4300). 
In contrast to 7(4230), the charmonium-like state 7(4300) 


is a very broad structure and its resonance parameters are 
7 ( 4300 ) = 4294 + 11 MeV and 77 ( 4300 ) = 201 + 55 MeV. 
Similarly to the case of \}/{2S)n^Ji^, there exist two 

different solutions to the cross section for e^e~ hcJi^n^. 
In Fig. 6, the fitting results and the resonance contributions 
are presented. Comparing the two solutions, we conclude 
that these two resonances interfere constructively and destruc¬ 
tively for Sol. 1 and Sol. 2, respectively. The cross section for 
e^e~ XcQ<^ is simulated by one resonance 7(4230), whose 
fitting results are presented in Fig. 7. Since the experimental 
data for this process are not abundant and have large errors, 
we also use the phase space of Xc 0 (^ to fit the exper¬ 

imental data. The fitting curve is presented in Fig. 7 with the 
X^ of 49.8, while the x^ of fitting the data with a resonance 
is 3.1. In addition, we further check the possibility of fitting 
the experimental data with the sum of the phase space and the 
resonance; the fitted is almost the same as that of fitting the 
data only with a resonance and hence the contribution from 
the phase space is ignorable, which indicates that fitting the 
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cross sections of e^e" XcOOJ with a resonance is more rea¬ 
sonable. On the experimental side, we expect that additional 
measurements from BESIII and the forthcoming Bellell for 
this process will provide a further restriction on the properties 
of F(4230). 

In the combined fit, we include the contributions of F(4230) 
in all three of these hidden charm production channels. Thus 
we can compare the branching ratios of 7(4230) decaying 
into if/(2S)7T^7:~, and XcO^^- From our combined 

fit, we find the product branching ratio for the XcO^^ chan¬ 
nel is ^ XcO^) - 2.2 + 0.6 eV, which 

is consistent with the value, 2.7 + 0.4 + 0.5 eV, reported 
by the BESIII Collaboration [10]. We list the ratios of 
S(T(4230) ^ il/{2S)n+7T-) and ^(7(4230) ^ hcTT^n) to 
S(7(4230) ^ Xc 0 (^) in Table III. In Scheme I, the ratio of 
S(7(4230) ^ ifi(2S)n^7T-) to S(7(4230) ^ Xcoco) is esti¬ 
mated to be 0.56 ± 0.25 and 0.59 + 0.26 for Sol. A and Sol. 
B, respectively, in which interference between 7(4230) and 
7(4360) is destructive as shown in Fig. 5. If interference 
is constructive, the ratio ” is relatively as small as 

0.13 + 0.12 and 0.14 + O.o'for Sol. B and Sol. D, respec¬ 
tively. 

As for the ratio of S(7(4230) —> hcn^n^) to S(7(4230) ^ 
Xc 0 <^)^ it is obtained as 0.07 + 0.05 and 3.21 + 0.97 for Sol. 
1 and Sol. 2, in which the interference between 7(4230) and 
7(4300) is constructive and destructive, respectively. In ad¬ 
dition, if 7(4230) is the iA(45), the spins of charm and anti¬ 
charm quarks in i/'(45) are parallel, while in he they are an¬ 
tiparallel. Thus, the il/{4-S) —> hcTi^n'' decay is a spin flip 
process, which should be suppressed in heavy quark effective 
theory. Thus, the ratio for Sol. 1, i.e., 0.07 + 0.05, is 

more favored as a physical solution compared with the one for 
Sol. 2. 


B. Scheme II 



(GeV) \/5 (GeV) 

FIG. 8: (color online). The same as Fig. 5 but in scheme II. 

In Scheme II, the cross section for hcn^ji" is 



FIG. 9: (color online). The same as Fig. 6 but in scheme II. 



FIG. 10: (color online). The same as Fig. 7 but in scheme II. 


simulated by a coherent sum of a Breit-Wigner function and 
2 —> 3 phase space as shown in Eq. (10), in which the 
cross section for goes up as the energy of 

the center-of-mass system is increased. The cross sections 
for —> \}/{2S)ji^n^ and e^e~ XcO^ are depicted by 

three Breit-Winger functions and one Breit-Winger function, 
respectively, which are the same as those in Scheme I. The 
fitted parameters are presented in Table IV, with which the 
;^^^/ndf is estimated to be 63.0/83 and is a bit larger than the 
one in Scheme I. In Scheme II, the resonance parameters of 
the 7(4230) are fitted to be 

»*r(4230) = 4220 + 8 MeV, 

rr( 4230 ) =43 + 9 MeV. 

The center value of the 7(4230) mass fitted in Scheme II is 
smaller than the one obtained in Scheme I, while the cen¬ 
ter value of the total width in Scheme II is much larger than 
Scheme I. 

Sine 7(4360) and 7(4660) are not involved in the other 
two processes, the resonance parameters of these two states 
are determined as mY(iYi 60 ) - 4360 ± 7 MeV, 7^(4360) = 
68 + 14 MeV, otf( 4360 ) = 4664 + 27 MeV, and 7^(4360) = 
93+21 MeV, respectively, by fitting the cross section for 
e^e~ —> il/{2S)n^Ji^ in Scheme II, which is consistent with 
the parameters in Scheme I within errors. The four solu¬ 
tions a - d to the cross section for e^e~ i/'(25)7r^7r“ 
in Scheme II are presented in Fig. 8. One finds that the 
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TABLE IV: The same as Table II but in Scheme II. 


final State 

0(2S );r+;r 

h^n^n 

XcQl^ 

Sol. a 

Sol. b Sol. c Sol. d 



™y(4230) 

4220 ± 8 



r y(4230) 

43 + 9 



P;*23„)S(7(4230) ^ /) 0.2 ± 0.3 

0.7+ 1.1 0.8 ±0.6 0.3 ±0.3 

0.3 ±0.2 

2.8 ±0.9 

™y(4360) 

4360 ± 7 



7 y(4360) 

68+14 



ry(:360)®(5'(4360) ^ /) 7.5 ± 1.1 

7.3 ± 1.5 9.0+ 1.2 6.1 ±2.2 



02 2.7 ±1.9 

3.9 ±0.7 4.1 ±0.7 2.4+ 1.9 



™y(4660) 

4664 ± 27 



7 K(4660) 

93 + 21 



ry(:660)®(5'(4660) ^ /) 5.9 ± 2.4 

2.0 ±0.6 6.1 ± 1.0 1.9 ±0.9 



03 1.0 ±1.9 

4.8+ 1.1 2.8+ 1.0 3.1 ±2.0 



Co 


12279 + 551 


04 


5.6 ±0.4 


Y^ndf 

63.0/83 




signal of 7(4230) is inconspicuous in the cross section for 
—> il/(2S)7:^7r^, which is caused by the relative large total 
width of 7(4230) fitted in Scheme II. In this case, the product 
rf(l230)®(^("^230) ^ iff(2S)n^n-) is fitted to be 0.2 + 0.3, 
0.7 +1.1, 0.8 + 0.6 and 0.3 + 0.3 for the four solutions, re¬ 
spectively. All these fitted products 7**4230)^(^(4230) 
tlr{2S)n^7t~) have large errors, which indicates that the data 
of the cross section for e^e~ —> xlj{2S)n ^are in disfavor 
with the broad 7(4230). 

The fitting results and the resonance contributions to the 
cross section for —> hcir^n^ are presented in Fig. 9. 

The cross section goes up with the increase of the energy of 
the center-of-mass system above 4.4 GeV and behaves as the 
phase space of e^e~ —> hcn^n~. Further experimental mea¬ 
surements of this cross section, especially above 4.4 GeV, will 
show the explicit tendency of the cross section, which will 
provide a cracial test for both Schemes. The cross section 
for —> Xc 0 (^ in Scheme II is presented in Fig. 10. The 

7(4230) with a lower mass and wide width also provides a 
good description for this cross section due to the large error of 
the experimental data. 

As in Scheme I, the ratios of S(7(4230) —> \p{2S)n^7i'~) 
and S(7(4230) ^ hcn^n^) to S(7(4230) ^ 7'cow) can be 
estimated by performing the combined fit and are listed In Ta¬ 
ble III. The ratio ^ has a large error and is compatible 

with zero. The large error dominantly comes from the error of 
the branching ratio of 7(4230) —> i/'(25 )7r^7r“. As for the ratio 
^ it is fitted to be 0.11 + 0.04, which is consistent with 
the expectation of heavy quark effective theory when taking 
7(4230) as ^(AS). 


V. SUMMARY 

Being stimulated by the anomalous mass gaps of the S- 
wave charmonia and the similarity between the charmonium 
and bottomonium families, we have predicted a missing S - 
wave vector charmonium state, xjj{AS ), located near 4.2 GeV 
in Ref. [3]. Theoretical estimates based on the quark pair 
creation model indicate that this predicted ) is a narrow 
state [3]. Since on the experimental side a narrow structure 
near 4.2 GeV has been observed in the hcn^ji^ [9] 

and e^e~ Xc 0 (^ processes [1], we have attributed the struc¬ 
ture in both processes to the missing (A(45). The search for the 
signal of i/'(45) in other channels is an intriguing and urgent 
problem. 

Recently, more precise experimental measurements of the 
cross section for —> >f/{2S )ji^n^ have been performed by 

the Belle Collaboration [1]. Since other than the structures 
of 7(4360) and 7(4660), a number of bump events near 4.2 
GeV have been found in the same cross section, in this pa¬ 
per we have performed a fit to the cross section with three 
resonances to find one additional resonance, 7(4230), with 
m - 4243 MeV and F = 16 + 31 MeV. The resonance 
parameters of the 7(4230) are consistent with our expecta¬ 
tion for the missing The upper limit of the prod¬ 

uct F**4230)®(^(4230) ^ \f/{2S)n^n~) is fitted to be 1.9 eV. 
With the assumption of 7(4230) as the ipiAS) and the theo¬ 
retical estimate of in a screening potential model, the 

branching ratio of i/^(45) —> \p{2S)n^n^ is evaluated to be 
S((A(45) ^ i/r(25)7r+;r-) < 3 x lO^l 

One should notice that this mass of the i/'(45) is above the 
threshold of a charmed or charmed-strange meson pair, which 
indicates that the dominant decay modes of xjjiAS) are open- 
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charm decays. As for the hidden charm decay processes, like 
i/^(45) —» i//(2S)7r^n^ and XcO(^, they can occur via charmed 
or charm-stranged meson loops. We have calculated the me¬ 
son loop contributions to the decay if/iAS) —> \p{2S)n^n'' 
to find that the branching ratio resulting from the meson 
loops is dependent on the parameters involved in calculations, 
but the branching ratio from our fit to the cross section for 
—» 1^(25 )7r^7r“ is understandable in a reasonable param¬ 
eter range, which also indicates that the 7(4230) observed in 
the —» ^{2S )n^n^ process could be a good candidate of 

the ij/(AS). 

To further test our conjecture that the structures near 4.2 
GeV in —> hcn^n^, —> il/{2S)n^n^, and 
XcO<JJ come from the same source, i.e., ip{A-S), we have per¬ 
formed a combined fit to these three hidden charm produc¬ 
tion processes in two different schemes. The mass and width 
of 1^(45) are fitting to be m^( 4 s) — 4234 + 5 MeV, r^( 4 s) = 
29 ±14 MeV and ~ 4220 ± 8 MeV, 7^(4^) = 43 ± 9 

MeV for Scheme I and Scheme II, respectively, as shown 
in Eqs. (II) and (12). Besides the resonance parameters 
of i/^(45), the products of the branching ratios of ifr(4S) 
^{2S)n^n^, hcTT^TT^, Xc 0 <^ the dilepton decay width of 
the i//(4-S ) are also determined by fitting the experimental data. 
The ratios of the branching ratios for the hidden dipion de¬ 
cay to the Xc 0 <^> decay mode are also estimated, which can be 
tested by further experimental measurements at BESIII and 


the forthcoming Bellell. 

Before closing this section, we would like to emphasize that 
the precise experimental measurements of the hidden-charm 
decay have provided some evidences of the existence of the 
missing (A(45) at the present moment. The future measure¬ 
ment for the hidden-charm decay will be able to further test 
our estimate of resonance parameters of iJ/iA-S ), in particular, 
the measurements of the cross section for —» hcn^n^ 

above 4.4 GeV can provide a further restriction on the prop¬ 
erties of >f/{4S). In addition, the precise measurements for the 
open-charm decay mode are also crucial for the properties of 
the missing ^{4S). We hope that precise measurements for 
both open-charm and hidden-charm decay modes will be car¬ 
ried out at BESIII and the forthcoming Bellell. 
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